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Abstract — An original multimode thru-reflection line (TRL) 
algorithm is used to derive the generalized scattering parameters 
of multimode two-port networks. Theoretical developments are 
detailed for calibration procedures based on two ports as well 
as multiple-port vector network analyzer (VNA) measurements. 
First-run experimental results demonstrate the validity of this 
technique. This method allows the experimental characterization 
of multiconductor transmission-line devices. It could also be 
used to characterize power coupling to undesired modes in 
monolithic microwave integrated circuit (MMIC) structures using 
conductor-backed coplanar waveguides. 

Index Terms — Calibration, multiconductor transmission lines, 
multimode waveguides, scattering parameters, scattering param- 
eters measurement. 



I. Introduction 

SCATTERING-PARAMETER measurements have become 
common and essential for microwave engineers. Vector 
network analyzers (VNA's) associated with efficient calibra- 
tion techniques such as the thru-reflection line (TRL) technique 
[1] offer almost rigorous measurements up to millimeter- wave 
range. However, conventional calibration techniques are based 
on a major assumption. Waveguides must propagate only one 
mode at the measurement reference planes. Most waveguides 
are designed so that only one dominant mode is present. If 
some evanescent modes are excited at discontinuities, refer- 
ence planes are chosen suitably far enough away so that their 
amplitude can be neglected. 

Nevertheless, multiconductor lines used in logic bus [2] or in 
monolithic microwave integrated circuit (MMIC) technology 
[3] are designed to drive several signals at the same time. As 
these lines propagate multiple modes, conventional calibration 
techniques generally fail. Similar situations are encountered 
with conductor-backed coplanar waveguides which propagate 
three dominant modes. One expects that only the symmetric 
coplanar mode will be excited in such structures. For this 
reason, coplanar probes are used and transitions are optimized 
in order to minimize undesired modes. However, it is well 
known from engineers that these type of circuits sometimes 
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exhibit strange behavior as compared to computer-aided design 
(CAD) tool predictions. This can arguably be attributed to 
microstrip-like mode excitation. 

In such situations, we need new calibration techniques able 
to handle multimode waveguides [4]. For conductor-backed 
coplanar lines, corresponding characterizations would at least 
allow the optimization of technologies and devices in limiting 
the influence of undesirable modes. For multiconductor lines, 
these new techniques could be used to build component 
libraries for interconnections, bends, crossing lines, via holes, 
and other discontinuities. 

With this in mind, we propose a generalization of the TRL 
calibration technique to multimode waveguides. In the first 
part of this paper, we will introduce the basic principle of this 
approach, which is based on a multiport scattering-parameter 
measurement. All measurements will be represented by a 
simple model. Generalized scattering and transfer parameters 
[5] which are used throughout this paper will be presented. 
The original multimode TRL algorithm will then be discussed. 
The theoretical developments will use notation similar to that 
used in the previous work of Eul and Schiek [1], which may 
be familiar to engineers involved in microwave calibration 
techniques. Finally, an experimental verification will be de- 
scribed in order to point out the validity of the multimode 
TRL approach. 

II. Theory 

A. Basic Principle of Multimode Calibration 

Multimode calibration intends to derive the generalized 
scattering parameters of multimode waveguide devices. As an 
example, we will consider the multiple conductor structure 
device under test (DUT), shown in Fig. 1. In general, this 
circuit may be a logic bus, coupled microstrip lines, or any 
coplanar waveguide. It may include crossing, open-end, bend, 
and other discontinuities, as well as active devices. We define 
two reference planes for which we assume that only N 
modes propagate. In order to extract the generalized scattering 
parameters, we need 27V* measurement channels, N on each 
of the two ports. With this in mind, the multimode DUT is 
connected to a VNA, as indicated in Fig. 2. N single-mode 
lines are used on each generalized port. For N equal to one, 
this simply describes the conventional single-mode calibration 



0018-9480/98510.00 © 1998 IEEE 



5- 1 k J 

SEGUINOT el at.: MULTIMODE TRL — A NEW CONCEPT IN MICROWAVE MEASUREMENTS 537 



Reference planes 




'Open end Crossover 

i PUT j 



I modes (N signal lines + ground)— 
Fig. L A multiple-conductor waveguide DUT. 



generalized 
portl 



Reference planes 



DUT 



generalized 
port 2 



N physical 

lines 
(N modes) 



1 multimode 

port 
(N modes) 



1 multimode 

port 
(N modes) 



N physical 

lines 
(N modes) 



B" 1 



Fig. 2. The basic setup for multimode TRL measurements. 

technique. Transitions A and B are quite unusual. They relate 
one multimode waveguide supporting N modes to N single- 
mode lines. A possible realization of transitions A and B 
will be presented in the experimental section. As indicated 
in Fig. 2, such measurements can be represented by simply 
cascading transfer matrices. Due to multimode propagation, 
we will require a generalized transfer and scattering matrices 
concept. 

B. Generalized Transfer and Scattering Matrices 

The DUT has two physical ports (1 and 2). Each port has N 
propagating modes. Such devices are commonly represented 
as 2iV-ports, each port being associated with one mode (see 
Fig. 3). Ports denoted li are located on physical port 1. More 
exactly, port li refers to the zth mode of generalized port 1. 
All incident (dji) and reflected (bji) modal-wave amplitudes 
are combined to define the generalized wave vectors Aj and 
Bj as follows: 



[Ai] = 







'hi" 


a J,2 


[By) = 








-a j, N - 




.bj, n . 



(1) 



Throughout this paper, uppercase letters refer to matrices, 
while corresponding lowercase letters refer to simple scalar 
matrix coefficients. Generalized waves are used to define the 
generalized transfer matrices [5] 

[2Mi]-S £][£]• « 

In this expression, Tij are (N x N) square-block submatrices 
which correspond to the DUT generalized transfer coefficients. 
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Fig. 3. Representations of a two-port multimode DUT. 

In the same way, we define the generalized scattering matrix 
(5) and coefficients (Sij) with 
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(3) 
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Using simple linear algebra, 5 and T matrices are related by 
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(4) 



This notation is close to the one used in [I] despite the 
large number of modes taken into account. In the single-mode 
case (N = 1), all block matrices (generalized coefficients) 
become single scalar so that these equations simplify to the 
conventional ones. 

As we have built a 2iV-port network, we require mea- 
surements of 27V-port scattering parameters. At this step, we 
propose two different approaches. Under some reasonable ap- 
proximations, a calibrated two-port VNA can be used to derive 
the required 2iV-port scattering parameters. This requires that 
all ports not connected to the VNA are loaded with perfectly 
matched terminations. Imperfect loads must be accounted for 
using an appropriate deembedding algorithm [6]. Moreover, 
calibration procedures are available to extract multiport 5- 
parameters from multiport VNA measurements [7]. Using 
these techniques, each measured transfer matrix M is related 
to the DUT T matrix by 



M = ATB 



-i 



(5) 



where A and J5" 1 are the transfer matrices of transitions A 
and By as described in Fig. 2. In this approach, the two-port 
VNA is first calibrated. The multimode TRL algorithm is then 
applied. As the measurement system is calibrated twice, this 
will be referred to as the "two-step calibration procedure." 

Using a multiport VNA (with 2N ports) could greatly 
ease the multimode TRL implementation in getting a "one- 
step calibration procedure." This is investigated below for the 
completeness although multiport VNA's are not commonly 
available. This paper has been organized so that readers who 
are not interested in the "one-step calibration procedure" can 
skip Section II-C. 
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Reference planes ' 
Fig. 4. A general model for multimode TRL measurements. 

C. A General Model for Multimode TRL Measurements 

The simple model proposed in Fig. 2 and (5) is based on the 
assumption that at some measurement planes "exact" multiport 
parameters are available. In this section, we will show the 
ability of such a simple model to describe all measurement 
situations involving a multiport VNA corresponding to Fig. 4. 
This model is a generalization of the single-mode double 
reflectometer model presented in [1]. The multiport VNA 
source and loads are successively switched to each port of 
the VNA. 4N measurements channels are used. Transitions 
A and B represent the linear mapping between measured 
complex-wave amplitudes (m^) and reference plane wave 
amplitudes (A{, Bi) as well as waves incoming from the 
generator and loads. Quantity mu (respectively, m2%) may be, 
for example, reflected (respectively, incident) complex-wave 
amplitudes measured on each single-mode input line (port 1 
in Fig. 2). We define measurements column vectors by 



m i,2 



m\y N ], » = !,.. -,4 (6) 



where the superscript (k) refers to one possible configuration 
of the switch. As in [1], it can be shown that measured values 
are only dependent on reference plane waves 



Meas^ 
Mea4 fc) 



= A 



(*) 



Mea4 fc) 
Measi fc) 



= B 



4 k) 

2?< fc > 



(7) 



where A and B are unknown (2N x 27V) matrices accounting 
for all imperfections of the VNA such as imperfect load, 
couplers, etc. Transitions from the VNA's ports and the 
reference planes are also taken into account by (7). The only 
assumptions which have been made concern the linearity of 
detectors and the repeatability of connections [1]. According 
to this model, each measurement must satisfy the equation 



(8) 



where T is the (2N x 2N) generalized transfer matrix of the 
DUT. Each configuration of the switch [referenced by the 
superscript (k)] gives one set of measurements. Combining 
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2N sets of measurements gives a (2N x 2N) matrix M, which 
will be called the measurements matrix 

Measi = [Uess[ l) Meas[ 2) Meaaj 2 ^]. 

i = 1 to N 



Measi 
Meas2 



Meas3 
Meas4 



-l 



(9) 



Using (7) and (8), matrix M is simply related to A and B by 
the same equation as (5). Thus, the same TRL algorithm will 
be developed for "one-step" as well as "two-step" calibration 
procedures. 

D. Multimode TRL Algorithm 

Starting from (5), each DUT of generalized transfer matrix 
N x will lead to the following measurements matrix M x \ 

M x = AN x B~ l N x = A~ l M x B. (10) 



We start calibration by measuring thru line (T) and line 
(L) standards. This choice is justified by the fact that only 
these components parameters are fully or partly known. This 
is especially true for multimode waveguides for which simple 
assumptions seem unrealistic. A line may be used instead of 
the thru line, as indicated by Eul and Schiek [1]. In our case, 
transfer parameters of thru (N\) and line (7V 2 ) standards are 
given by 



Ni = 



01 
[0 I 

N 2 =diag[exp(-7 1 Z),- 
exp(+7iZ), • 



•,exp(-7 N Z), 
•,exp(+7Ar/),] 



(11) 



where / (respectively, 0) denotes an (N x N) identity (respec- 
tively, null) matrix. 7* is the ith-mode complex propagation 
constant and / is the line length. As in [1], we define two 
similar matrices P and Q as follows: 

P = A- l QA. 

(12) 



P = N 2 N-i and Q == M 2 Mf 1 with 



We use the fact that similar matrices have identical eigen- 
values to define an eigenvalue matrix A as follows: 

eig(P) = eig(Q) = Xi A = diag(Ai), i = 1, N. (13) 

Since P is equal to #2* and A have identical eigenval- 
ues. This is used to derive the unknown propagation constants 
of the line 



exp(-7ii) = Xi. 



(14) 



Each mode is associated with two opposite propagation 
constants and a pair of inverse eigenvalues. Solving (14) may 
lead to numerical difficulties, as in the conventional TRL. 
The phase shift introduced by the line must be different from 
multiples of one half- wavelength. However, as this criteria can 
be fulfilled for one mode at one given frequency, it can fail for 
other modes with different propagation constants. We hope that 
multiple line standards can be used [8] in order to overcome 
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such problems and to improve the bandwidth and accuracy 
of the calibration, but this will not be investigated here. It 
is important to note that all propagation constants which are 
calculated by (14) must be correctly ordered in A so that N 
first (respectively, last) eigenvalues are related to N incident 
(reflected) waves. 

In fact, the algorithm which computes the 2N eigenvalues 
(13) cannot attribute a pair of eigenvalues to their specific 
existent modes (e.g., even and odd modes). In our algorithm, 
we have proceeded as follows. First, we identify N pairs of 
opposite eigenvalues. Each pair corresponds to one incident 
and the corresponding reflected mode. Then, the N eigenval- 
ues attributed to incident modes are sorted in ascending order. 
Therefore, the "first" mode (mode one) appears as the highest 
velocity mode. If the frequency band is limited, the phase shift 
does not exceed 180° so that the algorithm does not need any 
prior estimate of the mode's propagation constants. For highest 
frequencies, an estimate is required and can be provided by 
interpolating lower frequencies values. 

Up to this step, our notation is very close to that used by 
Eul and Schiek [1]. The derivation of propagation constants 
has been extended to multimode waveguide structures. As we 
are only concerned with thru-line-based calibration, we will 
now solve the problem in a different way from [1]. We use 
the following relation between P, Q, and A: 



A = X~ l PX = Y~ l QY 



(15) 



where the columns of X (respectively, Y) are composed of 
the eigenvectors of P (respectively, Q). Since eigenvectors 
are known except for an arbitrary constant, we can write X 
and Y as 



X=X Q f), 
Y=Y 0 6, 



0 = diag(ft) 
6 = diag(^), 



i = l,N. (16) 



Since P and Q are fully known, Xq and Y 0 can be computed 
from (15) while /? and S remain, for the moment, unknown. 
As P is diagonal, one possible solution for Xo is a unit matrix 
(Xq = J). Taking into account the relation between XY and A 



P = XY- l QYX~ l = A~ X QA. 
We use (16) to partly determine A as follows: 
A = YX~ l = YoKXq 1 = AoK 

A-=diag(AA) = 



Ki 0 
0 K 2 



(17) 



(18) 



where A$ is fully known while K remains unknown. The 
measurement of the thru is also used to partly derive B as 
follows: 



B = M^AN X = Aff 1 AoK = B 0 K< 



(19) 



At this step, the T and L standard measurements yield 
a partial determination of A and B while 2N unknowns 
in K remain to be solved. Any additional measurement of 
lines would lead to the same result. We need additional 
standard measurements to proceed with the calibration. We 
use a nontransmitting third standard (R). This one-port is 



characterized by its unknown generalized reflection coefficient 
denoted T (N x N matrix), which is defined in the same 
way as the generalized reflection coefficient S\i in (3). This 
generalized one-port (R), which is connected at both reference 
planes, is related to reference plane waves by 



Bi=rAi B 2 =TA 2 . 



(20) 



Corresponding measurements on generalized port 1 and 2 
(Fig. 2) yield two measured reflection coefficients denoted Ti 
and T 2 defined and related to T by 

Measi =TiMeas2 

Meas 4 — F 2 Meas3 (21a) 
r a = (A u r + A l2 ){A 2l T + A 22 y l 
T 2 = (B 22 T + B 2l )(B l2 T + Bn)" 1 . (21b) 

Equation (21a) is used only if a **one-step calibration 
procedure" is performed. In this case, the Measi are defined 
similarly as in (9), except that we only need N sets of 
measurements (instead of 2N) to derive a reflection coefficient. 
Using (18) and (19), it can readily be shown that (21b) can 
be rewritten as 



r = K 2 ^2^1 = 1 T\K 2 



-if 



(22) 



where the nonnormalized reflection coefficients f i and T 2 are 
known and defined by 

=[A) ( ii -riAo f 2i]~ 1 [Tii4o,22 - A), 12] 

?2 =[£o,22 ~ r , 2-Bo i i2F 1 |r , 2.Bo l ii ~ ^0,21 ]. (23) 

In these expressions Aq,^, B 0t ij denote (TV x N) subma- 
trices of Aq and Bo, as in (2) and (3). Solving the right side 
of (22) yields 

K 2 K^t x K 2 K^ =Lt l L- 1 =T 2 

WA« =f W (24) 

Where U denotes the ith term of diagonal (N x iV), matrix 
L and tk t ij is one coefficient of the matrix I*.. In (24), 
N 2 equations relate the N unknown diagonal terms of L. 
Nevertheless, solving (24), L is known except for one sign 
ambiguity (L = sLq with Lq fully known and 3 = ±1). This 
is used to relate K\ and K 2 as follows: 



K 2 = sLqKi , with 3 — ±1. 



(25) 



In the single-mode case, removing this sign ambiguity (3) 
is sufficient to complete the calibration. This is no longer 
true in the multimode case for which we have to provide an 
additional piece of information to solve (22). Any additional 
measurement of any arbitrary generalized one- or two-port 
network will not give additional information. Thus, the only 
way to proceed is to use the reciprocity of the reflection 
coefficient [3], [5], [9] which now has to be assumed for the 
third standard 



r = r t . 



(26) 
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Combining (22), (25), and (26), we obtain 



(27) 



It then appears that ifi is determined except for one 
arbitrary complex constant a and an (N - l)-order sign 
ambiguity 



10* 



(28) 



Inserting (28) and (25) in (22), we now get two estimates 
of the previously unknown reflection coefficient T valid in the 
multimode as well as in the single-mode cases 

r' =(aKior l ti(8aLQK lo ) = sK^TxLqK^ 

T" = (saL Q K l0 )- l t 2 (aK 10 ) = «(WTio)- l f 2*io- (29) 

We finally find as in the single-mode case (N = 1) [1] 
that an JVth-order sign ambiguity remains unsolved by the 
TRL algorithm. Comparing T' and T" in (29) does not help 
in removing the sign ambiguity. Thus, the user must provide 
an estimate of the R standard to proceed. Assuming that any 
sign ambiguity is removed, K is known except for the one- 
dimensional ambiguity (K = aKo). Inserting this in (10), 
we complete the calibration procedure and derive the transfer 
parameters of any DUT as follows: 

N x = a- l Ko l A^ l M x AoK 0 a = K^Aq l M x AoK 0 . (30) 

This is a generalization of (23) in [1]. Thus, a appears 
(as in previous works) as the one-dimensional ambiguity 
which is not removed by the TRL algorithm. However, this 
is sufficient for the determination of the DUT parameters. 
It can be shown that scattering parameters obtained using 
(4) and (30) are reciprocity normalized [5], [9]. In fact, 
reciprocity was essential to complete the calibration. At this 
step, important observations must be made. Equations (24) and 
(27) leading to the determination of K\ and K<i are redundant 
(except in the single-mode case). Thus, these equations may 
be solved using a least-square approximation. Furthermore, if 
unrepeatable connections of standards to the VNA ports leads 
to measurement uncertainties, V and T" will exhibit different 
values. This is used to define the following figure of merit F: 



\\T*-r"\\=F 



(3D 



where double vertical bars indicate any norm of a square 
matrix. This figure-of-merit must be as small as possible. This 
gives qualitative information about calibration precision which 
is evaluated without need of prior knowledge of the third 
standard scattering parameters. One may also observe that if 
the third standard is matched (T = 0), (22) becomes trivial and 
Ki cannot be determined. This might seem obvious. Since the 
first two standards are perfectly matched, the third one must 
be reflecting and preferably highly reflecting. Furthermore, 
solving (27) requires that at least one off-diagonal term of 
each line of T must be nonnull. Thus, the R standard must 
be designed to exhibit coupling between modes (I\j ^ 0 for 
i ^ j). For this reason, an "ideal" open or short cannot be used. 
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Fig. 5. Multimode TRL standards used for experiments. 
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Realizations other than a nontransmitting third standard (R) 
are possible. A transmitting network (N) would lead to a 
multimode TLN [1] calibration procedure. It can be shown 
that this would lead to conclusions similar to the generalized 
TRL procedure described above, and for this reason will 
not be presented here. Nevertheless, we must point out that 
TLN calibration offers the advantage of increased equation 
redundancy as compared to the TRL. 

At this step, we have demonstrated the theoretical feasibility 
of TRL calibration under the multimode propagation condition. 
We will now present our first experimental verifications. 

III. Experiments 

Experimental verification of the multimode TRL method has 
been made for two asymmetrically coupled microstrip lines. 
Because two dominant modes are propagated, two single- 
mode connections are used on each port, as indicated in 
Fig. 5. A 1 -cm-long line (L) standard which is approximately 
half-a-wavelength long at 5 GHz was used to allow TRL 
calibration up to 4 GHz. The third standard was realized 
by terminating and connecting together the two microstrip 
lines in the reference plane. Measurements were carried out at 
IEMN using a conventional two-port VNA. The VNA was first 
calibrated using the built-in "full two-port" menu. Connection 
of multimode standards to the VNA was made using microstrip 
to coaxial test fixtures. No corrections have been made to 
account for imperfect terminations used to measure 

the four-port scattering parameters. All results presented here 
correspond to first-run experiments. Our goal here is to prove 
the feasibility of multimode TRL and to illustrate potential 
experimental difficulties. All results have been derived without 
using any curve smoothing or averaging. 

Using our TRL algorithm, we have experimentally deter- 
mined the two propagation constants of the coupled microstrip 
lines. Mode "one" is chosen to have the lowest propagation 
constant. Slowing factors are presented in Fig. 6. Experimental 
values are compared to theoretical ones (dotted lines) obtained 
with the quasi-TEM model implemented in MDS. 1 The relative 
error between experimental and theoretical results is approxi- 
mately 4%. Thus, the accuracy of these results is obviously not 
sufficient to derive precise calibration. This may be attributed 
to poor repeatability of microstrip-to-coaxia! test fixtures as 
well as imperfect matching of 50-H loads. Better results can 

'"HP 851 SOB microwave and RF design systems, components catalog," 
Microwave Library Component, vol. 3, pp. 32-33, ch. 12. 
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Fig. 6. Slowing factors of coupled microstrip lines. Comparison of ex- 
perimental values derived from multimode TRL calibration to theoretical 
simulation. 
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Fig. 7. Generalized scattering parameters of a delay line, (a) Magnitude of 
the 16 S-parameters. (b) Phase of the four nonnull parameters appearing in 
(32). 



be expected using microwaves probes and accounting for 
imperfect terminations [7]. 

Completing the calibration, we have also derived the scat- 
tering parameters of a 1.5-cm-long delay line (see Fig. 7). 
Using our notation, parameter £21(1,1) * s > f° r example, the 
transmission coefficient of the lowest velocity mode (mode 
one) from port 1 to port 2. For an ideal transmission line, all 
parameters are null except the four parameters indicating the 
traveling of waves from one port to the other: 

S ij{k> ao = exp(- 7 *0, for * 5* 3 and k = 1, 2. (32) 

Our results indicate that the four above-mentioned parame- 
ters magnitudes are close to 0 dB, except around 5 GHz. The 
12 other parameters have negligeable amplitude (< -20 dB) 



up to 5 GHz and above 6 GHz. Measurement errors occurring 
between 5-6 GHz are inherent to TRL calibration. In this 
frequency range, the length of the line standard used in the 
calibration (1 cm) is close to a half-wavelength for the two 
propagating modes. This phenomenon which is consistent with 
conventional single-mode TRL is also observed in the phase 
of the scattering parameters. Measured phase shift of the two 
propagating modes [four scattering parameters of (32)] are in 
good agreement with the quasi-TEM model, 1 except in the 
5-6-GHz range. 

We have pointed out previously that off-diagonal terms 
of the reflecting standard must be nonnull. Therefore, we 
performed a full-wave analysis to determine the S-parameters 
of our R standard. We found a coupling between modes (I"^ 
and T2i) which was close to -20 dB. Such small reflection 
coefficients are sensitive to measurement uncertainties. At the 
same time, F12 and £21 (in 27) are in the same order of 
magnitude as Ti2 and ^1, and are used to solve (22). Such a 
situation might obviously be a source of calibration errors. 

Although the precision of these measurements is not suffi- 
cient, they do demonstrate the validity of our TRL algorithm. 
At the same time, these first-run results indicate that further 
efforts must be made to improve the accuracy of multimode 
TRL calibration. 



IV. CONCLUSION 

An original multimode TRL calibration procedure has been 
proposed. This procedure is based on measurements of scatter- 
ing parameters of a multiport circuit, including the multimode 
DUT. The corresponding algorithm is a generalization of 
the conventional single-mode TRL. The three multimode re- 
ciprocal standards used during calibration are a thru line, 
line, and reflect. The multimode TRL procedure leads to the 
determination of reciprocity normalized multimode scattering 
parameters. Using this original technique, we have experi- 
mentally determined propagation constants and generalized 
scattering parameters of asymmetric coupled microstrip lines. 
These first-run results exhibit insufficient accuracy for precise 
calibration; however, they demonstrate the validity of our 
algorithm. 

For MMIC's, multimode TRL could help us in analyzing the 
amount of energy coupled to undesired modes. For example, 
conductor-backed coplanar structures which propagate several 
dominant modes could be optimized to avoid any energy 
conversion from the coplanar mode to the microstrip one. For 
multiconductor transmission lines and logic buses, multimode 
TRL would allow precise characterization of crosstalk and 
discontinuities such as via-hole and crossing strips. With this 
in mind, we now have to improve this new method and extend 
it to other practical structures. This can be accomplished 
by taking into account multiple standards or symmetry of 
structures. A first possible problem arising from an improper 
choice of the reflecting standard has been outlined. We must 
also design a practical calibration kit for microwave probe- 
based measurements. A sensitivity analysis may also be helpful 
to point out the origin of calibration errors and to determine 
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the limits of this new technique. We hope that multimode TRL 
will open new paths in microwave measurements. 
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